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We analyze the dynamics of a continuously observed, damped, microwave driven solid state charge 
qubit. The qubit consists of a single electron in a double well potential, coupled to an oscillating elec- 
tric field, and which is continuously observed by a nearby point contact electrometer. The microwave 
field induces transitions between the qubit eigenstates, which have a profound effect on the detector 
output current. We show that useful information about the qubit dynamics, such as dephasing and 
relaxation rates, and the Rabi frequency, can be extracted from the DC detector conductance and 
the detector output noise power spectrum. We also demonstrate that these phenomena can be used 
for single shot electron spin readout, for spin based quantum information processing. 

PACS numbers: 78.70.Gq, 42.50.Lc, 03.67.Lx, 63.20.Kr 



Recently, rapid experimental progress in mesoscopic 
physics has meant that it is now possible to confine, ma- 
nipulate, and measure small numbers of electrons in sin- 
gle or coupled quantum dots |3, S S B • The mesoscopic 
environment of such confined electron systems may con- 
sist of phonons, electrons, and other electromagnetic de- 
grees of freedom. Thus these experiments are particularly 
interesting as they allow the complex interaction between 
such confined systems and their environment to be stud- 
ied at the single electron level. Furthermore, in view of 
the potential applications of such systems in solid state 
quantum information processing 5] , understanding these 
interactions is important, as it allows qubit decoherence 
mechanisms to be studied and accurately characterised. 

Of particular interest are coupled quantum dot (CQD) 
qubit systems driven by oscillating electric fields. The 
presence of a driving field resonant to the qubit en- 
ergy splitting can drive transitions into the excited state. 
Continuous measurement of such systems can reveal im- 
portant spectroscopic information about the qubit, such 
as the qubit splitting, Rabi frequency, and decoherence 
rates. Earlier work has focussed on current transport 
through driven, open CQD systems 6, 7]. In a recent 
experiment, a driven qubit comprising a single electron 
in a closed CQD system was continuously observed via a 
nearby quantum point contact (QPC) detector 

In this Letter, we theoretically analyse this system (see 
Fig. P). We account for the coupling of the CQD both 
to the QPC and to a generic bosonic environment, which 
may comprise of phonons or other electromagnetic de- 
grees of freedom 8j. Both the detector and the environ- 
ment contribute to the qubit relaxation and dephasing 
rates. We also demonstrate how resonant driving phe- 
nomena can be used for single-shot readout of the elec- 
tron spin. Our results are also relevant to other driven 
qubit systems, e.g. superconducting charge qubits 

A number of authors have considered the continuous 
measurement of undriven charge qubit systems by a QPC 
detector 0, 0? • Continuous measurement of driven 




FIG. 1: Schematic of a CQD charge qubit measured by a QPC 
under (a) non-resonant and (b) resonant microwave driving. 



superconducting flux qubits has also been considered T^ . 
In what follows, we adopt the quantum trajectories de- 
scription of the measurement process 0, ll^ [l5[. and 
generalize results obtained in previous work jly, |l3 on 
the measurement of undriven charge qubits using a QPC 
at arbitrary bias voltage. We first derive a master equa- 
tion (ME) for the dynamics of the damped, driven charge 
qubit system. We use solutions of the ME to determine 
the DC conductance and current power spectra of the 
QPC detector, and show how various qubit parameters 
can be extracted from measurements of these quantities. 
Finally, we describe our technique for spin readout. 
The model system we consider is shown in Fig. ^ for 



which the total Hamiltonian is given by H 

-f^drive H~ -^^meas H~ -^^leads H~ -^sb 



i^env, where {h= 1) 
H,y, = (-Aa, - ea,)/2 = -^a^'^ /2, (1) 

HdrWe — ^0 COs[(^ — r\)i\ (cos (5 Cr^ + sin ^ (Ja^) , (2) 
i^meas = ^ + ^Cr^ + H.C., (3) 



Pleads = ^^5,/e<^5,/e<^5,/c + ^ , (4) 

k q 

Hsh = crz^K.i{bl + bi), Henv = ^ujB,iblbi. (5) 

i i 

Here, i^sys is the bare qubit Hamiltonian, in which 
(Tx = \l){r\ + \r){l\ and a, = \l){l\ - |r)(r|, where \l) 
(|r)) denotes an electron state localized on the left (right) 
dot, (j) = vA^T? is the qubit energy splitting, and 
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\g){g\ — |e)(e|. i^drive corresponds to the driving 



field with frequency ujq = (j) — r]^ which may couple to 
both the (Tx and (Jz qubit operators, as parameterized by 
^- Hmeas dcuotcs the qubit-dctcctor coupling, in terms of 
the dimensionless tunnelling parameters T = \/2TigsgDT 
and V = ^27rgsgDX^ ^ = ^/V^^^rgsgD- We have as- 
sumed that the tunnelling amplitudes, T and x, and the 
densities of lead modes, gs and go, are approximately 
independent of k and q over the energy range where tun- 
nelling is allowed, i^ieads is the free Hamiltonian of the 
source and drain leads, where asj {(^dj) is the annihi- 
lation operator for an electron in the jth source (drain) 
mode, i^sb and iYenv correspond to a standard spin-boson 
coupling to a generic bath of bosons 18], where bi is the 
annihilation for the ith boson mode. 

To simplify the Hamiltonian, we first transform to an 
interaction picture defined by Hq = — a;ocri^V2 + ^ieads + 
i^env7 and make our first rotating wave approximation 
(RWA) and neglect terms oscillating rapidly compared 
to Qo and the microwave detuning, rj. Under this trans- 
formation, the time dependence is transferred to the in- 
teraction term and the Hamiltonian becomes 

m) = -yf^-^<7i'^^+Mt)®Y{t)+Bimz{t), (6) 

where Q = l^o sin(6> — S). We have defined the sys- 
tem operators Ai{t) = ^^e-'^^*Pn and Bi{t) = 
^^e-'^^^Qn, with ujn = 0, ±cJo, Pi = T + ucosOa, 

(e) 



n{uj)] + 27rJ{—u)n{—u), where Q{x) = (x + \x\)/2 is the 
ramp function, eV is the source-drain bias across the 
detector, J{u;) = ^i^l^ioo — ooB.i) is the bath spectral 
density, and n{(jo) is the thermal equilibrium Bose occu- 
pation number. To proceed, we make a second RWA, 
where we neglect terms in Eq.© rotating at a rate uoq. 
This RWA is justified in the limit of weak coupling to 
the detector and environment, uoq ^ 5'y,z(^n) ^26] . We 
also assume that the driving field is sufficiently weak that 
^y,z(^n + ^nO ^ ^y,z(^n), i-c that the noise spectra 
are slowly varying over frequencies of order Vt' . Finally, 
in order to treat dephasing within our perturbative tech- 
nique, we require that lim^^o J{^) uo^ where s > 1, i.e. 
that the bath is ohmic or superohmic at low frequencies. 

These approximations allow us to derive a ME which is 
valid for arbitrary source-drain bias and arbitrary bath 
temperature. However, in this Letter, we restrict our 
attention to the low bias (eV < (j)) and low temperature 
(kT <C (j)) regime, which has been probed in a recent 
1^]. In this case, the master equation for the qubit in the 
original interaction picture is given by 
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h:^V[a''^)]pi{t) + TrV[\g){e\]pi{t) 



(e) 



where r„ 



P2 = Pi = 



_ pdet 



(8) 



"penv 



is the pure dephasing rate. 



-u ^iii 0\e){g\^ Qi = cosOal , and Q2 = 
)l = — single) I . We have also defined the opera- 



tors Y{t) = XEkJ^- 

= Ei /^i(e"''^^'^^^l + H.c.) which act on the QPC and 
bosonic environment degrees of freedom, respectively. 

To derive the ME for the dynamics of the qubit 
alone, we further transform to a frame defined by 



with r^^* = 2u^cos^0Sy{0) and T^"^ = 2 cos^ i95'z(0). 
= r^^* + r^"^, is the relaxation rate, with T^^* = 
'^'a'^Dq^s^k + H.c.) and 1^^ cosH Sy {(j)) and T^"" = cosHSz{(l)). We have also 



in which all the dynamics 



HI, - 

are contained in the qubit- QPC and qubit-environment 
interaction terms. Then Hi'{t) = + 

Br{t) where Ar{t) = Y^nW e'^^^-^^^' Pnn' and 
Bi'(t) = E^^, e"'^'^"+'^-^^Qnn', for some operators Pnn' 
and Qri 



defined V[A]p = ApA^ - {A^ Ap + pA^A)/2. 

DC conductance The DC current through the QPC 
is related to the steady state occupation probability of 
the dot nearest the QPC, to first order in z^, by / = /q — 
5I{l\p{oo)\l) where /q corresponds to the current when 
the electron is localized in state |r) and 51 = AvTeV 
0,^3- Using Eq. 0, we can compute the steady state 
occupation of the left well, and the scaled conductance 
(M = 1 + (/ - Iq)/5I = 1 - {l\p{oo)\l)) is given by 



and 



0, ±n' where = ^/WTW- In 



this picture, the qubit density matrix, p, satisfies 

Pi'{t) = -Tr5,D,B{ f dt'[Hr(t),[Hr(t'),R(t')]]}. 
Jo 



M 



(9) 



(7) 



F^^ + F^/2 is the total dephasing rate. M 



We now make a Born-Markov approximation, setting 
the lower limit of the integral to —00 and R{t^) = 
pr{t) ® ps ® PD ® PB, where ps, Pd, and ps are 
equilibrium density matrices for the source, drain, and 
bath degrees of freedom. At this point, it is conve- 
nient to introduce the asymmetric quantum noise power 
spectra [H for the QPC and bath, Sy{oj) 
j^^dte^^''Tr[Y{t)Y{{))ps®PD\ = 6(cj-eV) + 6(cj + eV) 
and Sz{uj) = J^^dte''^^TT[Z{t)Z{0)pB] = 27rJ(cj)[l + 



where F^^/ 

is plotted as a function of e in Fig. |2(a)| Note that 
Fig. |2(a)| is in excellent qualitative agreement with re- 
cent experimental observations 4, 9, 10]. The absence of 
multi-photon resonances in Fig. |2(a)| is a consequence of 
the first RWA made above. 

The resonances occur at e = y^cj^ — A'^. Useful spec- 
troscopic information may be extracted from these reso- 
nant peaks. If the driving frequency, cjq, is known, and 
the Rabi frequency, is known independently (e.g. from 
observations of the time dependence of the detector cur- 
rent, as discussed below) then F^ and T^f can both be 
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FIG. 3: (a) Schematic of spin measurement scheme, (b) 
Conductance curves for spin-up and spin-down configurations 
showing distinct resonance peaks. 



FIG. 2: (a) Variation in conductance (scaled between and 1) where 



versus dot bias. Parameters Fc^ = Fr = = 0.03a;o and A = 
0.2a;o were taken to be constant over the entire range of e. (b) 
The peak height /i as a function of peak width 5e. Increasing 
microwave power scans from left to right. Also shown is the 
relationship between 5e and the Rabi frequency Qo- (c) Power 
spectrum, S{ix)). Inner peaks: — 7r/20, outer peaks: — 
7r/4. The tall peaks have detector limited dephasing, Fc^ = 
F^^* = 2z/^e'\/cos^ ^, whilst the short peaks are dominated 
by phonon dephasing, Fc^ = 0.03/0. Other parameters are 
eV/(t) = 0.5, zy = 0.1, r^o/^ = 0.5. 



determined. When cjq ^ A and assuming that F^, F^^/ 
and Vt do not vary significantly across the peak, from 
Eq. we find Q? ^ F^F^/ ((^eVF^, - 1), where Se is 
the half-width-half-maximum for the peak. Therefore, 
plotting Q against Se allows both F^ and F^^/ to be de- 
termined. However, in the absence of time-resolved mea- 
surements, Q may be unknown, because the relationship 
between the input microwave power and the electric field 
coupling to the qubit may be unknown. In this case, 
F(^/ can still be extracted by plotting the peak height, 
/i, against Je, for different values of the incident power. 
Again assuming coq ^ A and that F^, T^f and Q do not 
vary significantly across the peak, h and Se are related 
by 1/2 — F^//2(5e^, which is independent of the (un- 
known) quantity Q. These results are shown in figure 
|2(b)| Note that, for sufficiently weak driving, the peak 
width directly gives F^^/, while for stronger driving, the 
peak width is proportional to Q. 

Power Spectrum Further spectroscopic information 
may be obtained from the power spectrum of the cur- 
rent through the QPC, which is given by 6'(cc;) = 
2j^^dTG{T)e-'^\ where G(r) = E[I{t ^ T)I{t)] - 
E[I{t-\-r)]E[I{t)]^ is the current autocorrelation function, 
and . .] denotes the classical expectation. SJuj) can be 
computed using the same formalism as in [13, [I3 • In the 
low-bias regime the scaled (symmetric) power spectrum 
S{u) = S{u)/So (where So = 2elDC = e^T^eV is the 
shot noise background) is closely approximated by 



5H = 



Solo 



7q' 



70 



So 



2n^r^ + {2rf + 9?) Tr 

4^2pdet 

47?2i^2rdet(j^2p^2 ^ ^^^if + f}2)r,r^ + 4n2r^2) 



72, + (uj- ^'Y 72, + (a; + ^'Y 



(10) 



Thus the positions of the peaks yield the Rabi frequency 
Vt' ^ while the widths of the peaks contain other important 
spectroscopic information about the qubit. The height of 
the peaks at ±1^' is maximised when the external field 
is resonant with the qubit transition, r] = 0. Note that 
S{Q^) < 2 and thus the peak heights are no more than 
three times the shot noise background, as shown in Fig. 
|2(c)| The bound S{rt') = 2 can only be met when de- 
phasing is detector dominated, since extra dephasing due 
to the environment reduces the height of the peaks. 

Spin measurement Motivated by the preceding anal- 
ysis, we propose a method for single shot spin readout, 
using a microwave driving field and an inhomogeneous 
Zeeman splitting across the CQD. This splitting could 
be generated by an inhomogeneous magnetic field, or en- 
gineering different ^-factors in each dot. In this case, 
one spin configuration, say spin-down, may be made res- 
onant with the driving field, whilst the other, spin-up, 
is detuned by an amount 7^, and thus the spin can be 
determined by observing the current through the QPC 
[Fig. |3(a)| . Such a scheme is analogous to the method of 
spin readout via quantum jumps in atomic systems |2l| . 
We note that alternative methods for spin readout have 
also been proposed, via inhomogenous fields but without 
driving |22l. and by driving spin-flip transitions within a 
single dot [2^ . 

Spin can be included in the Hamiltonian by replacing 
(Jx and az in Eqs. (|1I5|) with generalized tunnelling and 
bias operators, i.e. ax ^l^^r} ^ ^u^ri + 
ni — rir where rii = c\^Ci^ + ^li^U denotes the number of 
electrons on site i. The inhomogenous Zeeman splitting 
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is included by adding the term 

^Zeeman = " ^ ^ BiS\^^ (11) 

to i^sys5 where Bi denotes the Zeeman sphtting on each 
site, and s\^^ = cJ^Ci^ — cJ^Ci;. We neglect any explicit 
coupling between spin up and spin down states (processes 
which can induce such transitions can be accounted for 
by a finite spin fiip time, which we discuss below). 

The Zeeman term amounts to a spin dependent bias 
between the dots, e^^ = e ± {Br — Bi)/2. The DC re- 
sponse of the detector is therefore shifted for each spin 
configuration, as shown in Fig. |3(b)| If Br — Bi > the 
resonant peaks are clearly resolved. In this case, if e is 
tuned such that the spin-down transition is resonant with 
the driving field, the currents for each spin configuration 
are approximately ^ Iq — and ^ Iq — 61 /2. If the 
detector is shot noise limited, the time taken to resolve 
these currents is rf/ = {1^^ - l]'^fl2e ^ 8P /m^. 
Thus = 4roi where tqi is the time taken to distin- 
guish the currents due to and 1 electrons on the dot 
adjacent to the QPC. It is believed that shot noise lim- 
ited QPC detectors will be available in the near future 
with Toi ~ 25 ns [2^, and so t^i ~ 100 ns should be pos- 
sible. This compares favourably with recently observed 
spin fiip times (Ti ^ 1 ms) in GaAs quantum dots 3]. 

The use of microwave driving for spin readout offers 
some advantages over other schemes where no driving is 
used. Firstly, a relatively small differential Zeeman split- 
ting is needed. In order to clearly resolve the resonant 
peaks, we require B^ — Bi> 5e. The lower bound on the 
peak width is given by 5e > F^^/. So for a charge dephas- 
ing rate of F^/ 10^ s"^ 2], we require Br - Bi > 0.07 
jneV. If one attempts readout via an inhomogenous Zee- 
man splitting but without driving [2||, then to obtain a 
comparable signal-to-noise, the differential Zeeman split- 
ting must be larger than the central transition region in 
Fig. [3(b)l i.e. Br - Bi > max(A, kT). For T = 100 mK, 
Br - Bi >9 fieV is required. Thus in GaAs {g = 0.44), 
with a uniform field of 1 T, our scheme requires a ^-factor 
variation between the dots of Ag/g ~ 0.3 %, whereas 
without driving, one would require Ag/g ~ 35 %. 

Secondly, the scheme can be used in such a way that 
a definite signal is always obtained for both spin con- 
figurations. Switching the microwave frequency first on 
resonance with the spin-down transition, and then on res- 
onance with the spin-up transition yields a definite signal 
indicating the spin state. That is, when the driving fre- 
quency is switched, the DC conductance will increase if 
the electron is spin-up, and will decrease if it is spin- 
down. This is in contrast to other measurement schemes 
in which a definite si gna l is only registered for one spin 
configuration, [23, |2J, I23 , with the other state indicated 
only by the lack of a signal. This has the advantage of 
eliminating false negative signals, where an "no-signal" 



event is erroneously recorded as evidence for a particu- 
lar spin configuration. Thus the resulting measurement 
fidelity should be improved. 

In summary, we have analysed the dynamics of a con- 
tinuously observed, driven solid state qubit, coupled to 
a generic bosonic environment. Both the environment 
and the coupling to the detector contribute to the qubit 
relaxation and dephasing rates. Useful spectroscopic in- 
formation, in particular the dephasing rate F^^/, can be 
extracted from DC measurements of the detector output 
current alone, even when the coupling between the mi- 
crowave field and the qubit is unknown. If the power 
spectrum of the detector output noise can also be mea- 
sured, then the relaxation rate F^ and Rabi frequency Q 
can also be determined. We have also proposed a sin- 
gle shot spin readout technique using microwave driving, 
which offers advantages over existing schemes and can be 
implemented with current technology. 
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